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We extend Mode-Coupling Theory (MCT) to inhomogeneous situations, relevant for supercooled 
liquid in pores, close to a surface, or in an external field. We compute the response of the dynamical 
structure factor to a static inhomogeneous external potential and provide the first direct evidence 
that the standard formulation of MCT is associated with a diverging length scale. We find in 
particular that the so called "cages" are in fact extended objects. Although close to the transition 
the dynamic length grows as \T — Tc|~'^''* in both the f3 and a regimes, our results suggest that the 
fractal dimension of correlated clusters is larger in the a regime. We also derive inhomogeneous 
MCT equations valid to second order in gradients. 
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It is becoming increasingly clear that the viscous slow- 
ing down of supercooled liquids, jammed colloids or gran- 
ular assemblies is accompanied by a growing dynamic 
length scale, whereas all static correlation functions re- 
main short-ranged. This somewhat unusual scenario, 
suggested by the experimental discovery of strong dy- 
namical heterogeneities in glass- formers jl| , ha_s_been_sub- 



stantiated by detailed numerical simulations [_,_,__,_,, 
explicit solution of simplified models 0, ^3 ^.nd very re- 
cent direct experiments 0, where four-point spatio- 
temporal correlators are measured. From a theoretical 
point of view, our understanding of supercooled liquids 
owes much to the Mode-Coupling Theory (MCT) of the 
glass transition. Although approximate in nature, MCT 
has achieved many qualitative and quantitative successes 
in explaining various experimental and numerical results 
[ill ll^ . In spite of early insights , the freezing pre- 
dicted by MCT was repeatedly argued to be a small 
scale caging phenomenon, without any diverging collec- 
tive length scale. This, however, is rather surprising from 
a physical point of view, since one expects on general 
grounds that a diverging relaxation time should involve 
an infinite number of particles [l5| . Building upon the 
important work of Franz and Parisi 0|, two of us (BB) 
[13 suggested a way to reconcile MCT with physical in- 
tuition. Within a field theory formulation of MCT, BB 
showed that the four-point density correlation function is 
given by the so-called 'ladder' diagrams that indeed lead, 
upon resummation, to a diverging dynamical correlation 
length and spatio-temporal scaling laws. BB also pro- 
posed a Ginzburg criterion that delineates the region of 
validity of MCT, which breaks down in low dimensions. 
Still, the field theory language used in [3| is not triv- 
ially related to the standard, liquid theory formulation 
of MCT . Indeed recent work has shown that the 
field theory is laden with subtleties ll^ , in par- 
ticular related to the Fluctuation-Dissipation relation. 
The aim of the present letter is twofold. First, we show 
how the results of BB may be recovered and extended 



to obtain testable, quantitative predictions on absolute 
dynamic length scales, entirely within via the standard, 
projection-operator based MCT Onr detailed anal- 
ysis predicts a remarkable scaling behavior that has im- 
plications for the geometry of dynamic heterogeneities. 
Second, our formulation generalizes MCT to spatially in- 
homogeneous situations. 

In order to proceed, we consider an atomic fluid sub- 
ject to an arbitrary external potential C/(x), such that 
the equilibrium averages (e.g. the static density /^(x)) 
vary in space. This opens a route to an MCT analy- 
sis of various situations of experimental interest, such 
as liquids in pores [U, close to a wall or a free 
surface, or sedimentation effects. The relationship with 
the results of BB will be obtained using this inhomoge- 
neous MCT formalism to compute the response of the 
dynamical structure factor to a localized external po- 
tential. In the limit where the wave vector associated 
with the external field tends to zero, a connection to 
the 4-point correlator of BB emerges. The dynami- 
cal quantities of interest are the density fluctuations Sp 
and the currents J, defined in Fourier space as dpq = 

Eti e*'' ""' - <^(q) and Jq = Eti q • P^e''!---/™, where 
(/i(q) = 'i- / N (^J2iLi ^^'^'^''^ ■ Following standard proce- 
dures based on the Mori-Zwanzig formalism 0, E^f . 
one can establish the following exact equation of mo- 
tion for the dynamic structure factor F{qi,q2;t) = 



J rfq'i^dt' M(qi,q;;t-i')^J^(q'i,q2;<') = 0, 



(1) 



where n^{qi,q[) = ^Slq, . q'^0(q, _ q'i)5-l(qi, qi) 
(with S ^(qi,q2) the inverse operator of i^(qi,q2,t — 
0)), and M(qi,q2;t) is a memory kernel which can be 
expressed in terms of the fluctuating part of the force. 
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Factorization approximations, analogous to those used 
in standard MCT, reduce this memory kernel to two- 
body correlation functions and yield an inhomogeneous 
mode-coupling theory (IMCT). The general IMCT equa- 
tions are rather cumbersome and will be presented else- 
where 12^. In the limit C/(x) 0, they reduce to 
the standard MCT equations. In order to obtain some- 
what tractable expressions, one can consider weakly in- 
homogeneous situations U (x) <C kBT, such that one 
can expand all quantities to first order in the pertur- 
bation. The aim is to compute the sensitivity of the 
dynamical structure factor to a small perturbation of ar- 
bitrary spatial structure, in particular localized perturba- 
tions, which wc can always decompose in Fourier modes: 

where Xqo(qi'^) « ^^'''^suTq!^)^''*^ ^u=o response of 

the dynamical structure factor to a static external per- 
turbation in Fourier space. For a perturbation localized 



at the origin, f/(x) = Uo5{x), one finds 6F{qi,y,t) = 
J dqos^'^°'^Xqo{^i^t)- This susceptibility is akin (al- 
though not exactly related) to a three-point density cor- 
relation function in the absence of the perturbation. Al- 
though quite different from the four-point functions con- 
sidered previously in the literature, Xqo(qi! ^) expected 
to reveal the existence of a dynamical correlation length 
of the homogeneous liquid (see Q for the particular case 
qo = 0). Indeed, XqelQii^) nieasures the influence of a 
density fluctuation at a given point in space on the dy- 
namics elsewhere. Within MCT, this three-point correla- 
tion turns out to have exactly the same critical behaviour 
as the one obtained by BB. As shown below, the deep un- 
derlying theoretical reason for such a coincidence is that a 
certain linear operator becomes critical at the transition 
(see [25! for a diagrammatic explanation). Differentiat- 
ing Eq. iQJ with respect to [/(qo), the final equation for 
the susceptibility XqQ(qi,t) reads: 



dk 



mS{qi) ^" ' Jo ' dt' Jo "^|qi+qo| 

-Kqi; k, qi - k)viq, + qo; qi - k, qo + k)xqo(k, t ~ t')^^o(|qi - k|, i - i')^MHl±S^ 



(2) 



dt' 



'^qolqi.O 



r 



where v is the usual MCT vertex i;(q;ki,k2) = q • 
kinc(/ci) -|- q • k2nc(fc2); Mq is the MCT memory 
kernel: Moiq,t) = knTp/2m J dk/{2'K)^v^{q;k,q - 
k)Fo(fc, t)Fo(|k — q|, t) [23. The source term S ^j, (q i , t), 
whose precise form will be presented elsewhere |24j . de- 
pends on Fo(k,t' < t) and on static four and five point 
density correlations; the value of the dynamic length scale 
and the critical properties of Xqo are however indepen- 
dent of the precise form of this source term. The above 
equation is of the type Cq„ Xqo — '^qo ' where is a cer- 
tain linear operator, the structure of which - in particu- 
lar its smallest eigenvalue - contains the information we 
want to extract. One should first note that in the limit 
qo = 0, the operator £0 simply encodes the change of 
the dynamic structure factor when the coupling constant 
(i.e. the density or the temperature) is shifted uniformly 
in space. This remark allows one to compute xo(k, t) 
from standard MCT results in the f3 and a regimes: 



Xo(k,t) 



7= — 9P % = U, — 



1 



-ga.k 



t 



-l/2a-l/26 



(3) 

where we use standard MCT notations [ill i^ik) is 
the structure factor and h(k) is the critical amplitude 
and e is the distance, in the liquid phase, between the 



coupling constant and its critical value). The behavior 
of the scaling functions at large and small argument can 
be found directly by analyzing Eq. Q or by scaling: 
in the early (3 regime u = t/rp the e dependence 
should drop out, hence g/s{0,u) oc u". The matching 
between a and /? regime implies 17/3 (0,u) oc at large 
u, whereas ga,k{u) oc S{k)h{k)u^ at small u. How is 
the behaviour modified for small modulations with qo ^ 
0? The derivation is simple for Xqo(qii°*2) the glass 
phase, where straightforward manipulations of Eq. ^ 
allows one to show that it satisfies the matrix equation 
(/ — M) ■ Xqo(qi'°o) ~ ■5qo where Sq^ is a source term 
that is regular and of order one in the limit Q and: 



pG(qi)G(qo + qi)S'(qi - q2) 



(2^)3gi|qi +qo| 



^^qo(qi'q2) 

/qi-q2''^(qi;q2,qi - q2)w(qi + qo;qi - q2,qo 

where G(qi) 



q2), 
(4) 

S'(qi)(l — /qj and /q is the non-ergodic 
parameter. Interestingly, the matrix M is exactly the 
same as the one obtained from the resummation of the 
ladder diagrams in the field theoretical framework of 
BB. Note also that the source Sq^ is irrelevant pro- 
vided it is not orthogonal to the lowest eigenvector of 
M . For go = and e < 0, the largest eigenvalue of 
M is, as shown by Gotze, A = 1 — 0(^/— e) |2^ and 
its right eigenvector is S{k)h{k). The correction 6\ to 



this eigenvalue at go ~* can be computed by pertur- 
bation theory. By symmetry, one expects that in gen- 
eral i5A = — r^Q, where F is a certain coefficient, lead- 
ing to Xqo(k,oo) ~ S{k)h{k)/{^\+Vql) In the 

schematic limit where 5'(q) is sharply peaked around 
q ~ Q, with a small width AQ, one can compute T ex- 
actly; one finds that F is positive and oc AQ~^. For a 
more realistic shape of S{q), e.g. hard sphere structure 
computed within the Percus-Yevick approximation at the 
mode-coupling density 0c = 0.515, we have determined 
5X numerically to extract F — 0.072(t~^. For more real- 
istic hard-sphere structure factors F may be as large as 
F = 0.3ct~^. We have not been able to show in full gen- 
erality that F should always be positive. A negative F 
would predict a remarkable 'modulated' glass transition, 
with the non-ergodic factor displaying periodic oscilla- 
tions in space |27l |. 

The analysis of the /3 and a regime in the liquid phase is 
more involved. As previously, the operator that becomes 
critical at the transition turns out to be the same as the 
one considered in BB (an d whose inversion leads to ladder 
diagrams). We find Hi that (/ - M^) ■ Xqo(^^^^^°) = 
Sqi^(te^/^°), where Mc is the matrix M at the transition 
and the new source term is of order one in the limit qo 
0. As a consequence one finds jl^: 

where 5^ satisfies qQgi3{qQ/ y/e.te^/'^"-) = (-y/e + 
F(7Q)(/|Sq^(te^/^")) with (/| is the left eigenvector con- 
jugated to S{k)h{k). The analysis of Xqo(k, i) in the 
a-relaxation regime turns out to be more subtle than an- 
ticipated in BB and will be detailed elsewhere We 
have established that for small go and fixed e, 

with S a certain regular fimction with S(0) 7^ and 
S(t; ^ 1) ^ 1/v such that Xqo behaves as q^^ for large 
go, independently of e. Also, ga,k{u ^ 1) = S{k)h{k)u'', 
as to match the /3 regime, and ga{u » l,k) 0. All 
those analytical results are in full agreement with the nu- 
merical solution of the schematic version of IMCT equa- 
tions (see Fig. 1). Note that the scaling variable is still 
gpE"^/^ in that regime, rather than go£~^ as surmised 
in BB . The physical consequence of the above anal- 
ysis is the existence of a unique diverging length scale 
^ ^ \/r|£|~^/^ that rules the response of the system to a 
space-inhomogeneous perturbation and hence of the spa- 
tial dynamic correlations. The analysis of the early /3 
regime where t ^ jel"^/^" shows that this length in 
fact first increases as i"/^ and then saturates at ^. Fur- 
thermore, Eq. 13) indicates that although the integrated 
dynamic correlation Xqo=oO^-i ^) increases in the a regime 

as £(.b-a)/2a-i.b (^^^^^ ^-1/2 f^^. ^ = ^-^ tO £"1 for t = T^) thc 
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FIG. 1: Numerical solution of the schematic IMCT equa- 
tions for T > Tc, where all k dependence is neglected (see 
for details). Main plot: Xqo(^) different qo as a 
function of time, and for e = 10~^. From top to bottom: 
go = 0, 0.06, 0.2, 0.4, 1.. Note that the shape of Xqo (i) in the a 
regime is independent of qo, as predicted by Eq. 0. We have 
in fact checked that the predicted scaling is very well obeyed 
in that region. Inset: exmax(qo) = eXqo(^ = t^) as a function 
of qoe~^^^ in log-log, for different go's and e's. Note the g^^* 
behaviour for large goe^^^*, as indicated by the dashed line. 



dynamic length scale itself remains fixed at ^. Interest- 
ingly, this suggests that while keeping a fixed extension 
^, the (fractal) geometrical structures c arry ing the dy- 
namic correlations significantly 'fatten' [28j between Tp 
(where the structures could correspond to the strings re- 
ported in recent simulations (S. i23]) and Tq, where more 
compact structures are expected, as indeed suggested by 
the results of [s^- For ^ i <C Tq., we expect a cross- 
over between dense and dilute structures at a new, time 
dependent crossover length . 

Starting from the general IMCT equation Eq. one 
could have chosen to follow a slightly different path and 
only assume that the length scale £ of the imposed inho- 
mogeneities is large, and perform a gradient expansion 
to order to obtain an equation on the space depen- 
dent structure factor, F(q, r, t). This space-dependent 
Ginzburg-Landau like MCT equation has one part that 
is identical to standard MCT equation (with space depen- 
dent coefficients) plus non linear contributions containing 
a V^-F term and, interestingly, a Burgers non-linear term 
of the form (VF)^ When inhomogeneities are small, 
F(q, r, t) = Fo(q, t) + x(q, r, t) with x < Fo, these equa- 
tions become identical to Eq. ^ above. In the schematic 
limit where all wave-vector dependencies drop off, our 
equation coincides with the gradient expansion obtained 
by Franz for the p-spin glass model in the Kac limit [sJl ■ 
In summary, we have shown how to extend the stan- 
dard framework of MCT to inhomogeneous situations. 
This allows us to compute the response of the dynam- 
ical structure factor to spatial perturbations. The case 



4 



of a localized perturbation at the origin is particularly 
interesting: it shows directly that the dynamical struc- 
ture factor is affected on a dynamic length scale ^ that 
diverges as \e\~'^ as the Mode- Coupling transition is ap- 
proached, with v = 1/4: and a prefactor that can be com- 
puted numerically. This length scale governs both the 
(3 and a relaxation regimes, showing that the standard 
interpretation of the /3-regime as the vibrations of par- 
ticles trapped in independent cages formed by nearest 
neighbors is somewhat misleading: as the MCT transi- 
tion is approached, ^ grows and the 'cages' become more 
and more collective. Our results suggest an interesting 
scenario where the geometrical structures responsible for 
dynamic fluctuations thicken with time. Note that ^, as 
in ordinary critical phenomena, only diverges at Tc, re- 
flecting the critical fragility of the system right at the 
transition. It is therefore clearly distinct from the di- 
verging viscous length y/rjTa which sets the scale below 
which the liquid sustains shear waves |32l |. which is in- 
finite in the whole glass phase. Our IMCT equations 
should be very useful to investigate inherently inhomo- 
geneous physical situations. Also, spatial critical fluctu- 
ations arc expected to play a major role in low dimen- 
sions, as for usual critical phenomena. These should lead 
to non trivial values of critical exponents (such as ly) 
and be involved in the breakdown of the Stokes-Einstein 
relation between viscosity and diffusion However, 
these critical fluctuations should also interfere with ac- 
tivated 'droplet' fluctuations that are expected to smear 
out the MCT transition in finite dimensions j^^l- The 
details of the way these two type of phenomena interact 
and lead to the observed crossover between an MCT-like 
regime and an activated regime in supercooled liquids is, 
in our opinion, one of the most crucial open theoretical 
questions. Mixing the present framework with the ex- 
tended Mode-Coupling scheme recently proposed in Q 
might be a promising path. Finally, we want to stress 
that generalized susceptibilities such as (k, t) offer a 
direct way to study correlations between dynamics and 
local structural fluctuations They provide a comple- 
mentary and more direct physical information than the 
4-point correlations studied in |^, |M |^ . We notice that 
Xqo (k, t) can be measured using state of the art molecu- 
lar or Brownian dynamics simulation. Experimentally, it 
could be accessible in colloids by use of an optical tweezer 
array imposing a periodic dielectric force on the particles. 

We warmly thank L. Bcrthier, S. Franz, W. Kob, V. 
Krakoviack, G. Tarjus, M. Wyart and F. Zamponi. GB 
is partially supported by EU contract HPRN-CT-2202- 
00307 (DYGLAGEMEM) and DRR by grants NSF CHE- 
0505939 and NSF DMR-0403997. 



[1] M. D. Ediger, Ann. Rev. Phys. Chem. 51, 99 (2000). 



[2] 
[3] 

[4; 

[5 



[9: 

[10 

[11 

[12; 

[is: 



[14 



[is: 
[16: 
[17: 

[is: 

[19 

[20 
[21 



[22: 
[23: 

[24: 

[25: 

[26' 
[27 

[28: 

[29 
[30 

[31 

[32' 
[33: 

[34] 



M. Hurley, P. Harrowell, Phys. Rev. E 52, 1694 (1995). 
R. Yamamoto and A. Onuki, Phys. Rev. Lett. 81, 4915 
(1998). 

G. Parisi, J. Phys. Chem. B. 103 4128 (1999). 

C. Bennemann, C. Donati, J. Bashnagel, S. C. Glotzer, 
Nature, 399, 246 (1999), S. C. Glotzer, J. Non-Cryst. 
Solids 274, 342 (2000), see also [3. 

L. Berthier, Phys. Rev. E 69, 020201 (2004); S. White- 
lam, L. Berthier, J. P. Garrahan, Phys. Rev. Lett. 92, 
185705 (2004). 

O. Dauchot, G. Marty and G. Biroli, Phys. Rev. Lett. 95 
265701 (2005). 

L. Berthier, G. Biroli, J. -P. Bouchaud, L. Cipelletti, D. 
El Masri, D. L 'Hot^ F. Ladieu, M. Pierno, Science 310, 
1797 (2005), and ^ 

F. Ritort, P. SoUich, Adv. Phys. 52, 219 (2003). 

J. P. Garrahan, D. Chandler, Phys. Rev. Lett. 89, 035704 
(2002). 

W. Gotze, L. Sjogren, Rep. Prog. Phys. 55 241 (1992). 
W. Gotze, Condensed Matter Physics, 1, 873 (1998). 
T.R. Kirkpatrick and D. Thirumalai, Phys. Rev. A 37, 
4439 (1988); T.R. Kirkpatrick, P.G. Wolynes, Phys. Rev. 
B 36, 8552 (1987) 

S. Franz, G. Parisi, J. Phys.: Condens. Matter 12, 6335 
(2000); C. Donati, S. Franz, G. Parisi, S. C. Glotzer, J. 
Non-Cryst. Sol., 307, 215 (2002). 
see A. Montanari, G. Semerjian, cond-mat/0603018 

G. Biroh, J. P. Bouchaud, Europhys. Lett. 67, 21 (2004). 
K. Miyazaki, D. R. Reichman, J. Phys. A: Math.Gen. 38, 

L343 (2005). 

A. Andreanov, G. Biroli, A. Lefevre, |cond-mat/0510669| 
M. E. Gates and S. Ramaswamy, Phys. Rev. Lett. 96, 
135701 (2006). 

G. Biroli, J. P. Bouchaud, in preparation. 

see V. Krakoviack, Phys. Rev. Lett. 94 065703 (2005) 

for MCT in pores, but for averaged quantities for which 

explicit inhomogeneities have disappeared. 

P. Scheidler, W. Kob, K. Binder, J. Phys. Chem. B, 108, 

6673 (2004) 

D. R. Reichman and P. Charbonneau, J. Stat. Mech. 
P05013 (2005). 

G. Biroli, J.-P. Bouchaud, K. Miyazaki, D. R. Reichman, 
in preparation. 

L. Berthier, G. Biroh, J.-P. Bouchaud, W. Kob, K. 

Miyazaki, D. R. Reichman, forthcoming. 

W. Gotze, Z. Phys. B 60, 195 (1985) 

This modulated glass bears some resemblance with, but 

is different from, the stripe glass of J. Schmalian, P. G. 

Wolynes, Phys. Rev. Lett. 85, 836 (2000). 

Our results suggest that in high dimensions, the fractal 

dimension of the dynamic clusters increases from df = 2 

in the (3 regime to d/ = 4 in the a regime. 

M. Vogel, B. Doliwa, A. Heuer, S.C. Glotzer, J. Chem. 

Phys. 120, 4404 (2004). 

G. A. Appignanesi, J. A. Rodriguez Fris, R. A. Montani, 

and W. Kob, Phys. Rev. Lett. 96, 057801 (2006) 

S. Franz and F. L. TonineUi, J. Phys. A: Math. Gen. 37 

(2004) 7433 and private communication. 

R. Ahluwalia, S. P. Das, Phys. Rev. E 57, 5771 (1998). 

T. R. Kirkpatrick, D. Thirumalai, P. G. Wolynes, Phys. 

Rev. A, 40 (1989) 1045; J. P. Bouchaud and G. Biroh, 

J. Chem. Phys. 121, 7347 (2004). 

P. Mayer, K. Miyazaki, D. R. Reichman, 
,cond-mat/0602248, 



